Effects of the neutral point of dust charge 
in plasma sheath 
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Abstract: We numerically investigate the dust charging in the sheath, by using the usual fluid 
approximation, it is extended to include self consistently the dust charge variation. The grain 
charge becomes a new self consistent dynamic variable, it is found that dust grains are first 
charged negatively at the sheath edge and then begin to be charged positively in the sheath. 
The numerical results show that the initial velocity of dust charged grains, sized grains and the 
density ratio of ion to electron at the sheath edge have affected the neutral point of dust charge. 
Moreover, the neutral point of dust charge affects the spatial distribution of dust density. 
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1. INTRODUCTION 

The study of the dusty plasma sheath is one of the im- 
portant problems in the experiment equipments of plasma. 
So the correlative works both of experiment and theory 
are developed widely during the past several years [1-10]. 
Recently, the sheath of dusty plasma has been studied with 
constant charge dust grains [3,4] . Many authors have inves- 
tigated the problems to include the effects of dust charge 
variation [1, 9-12], the charges of dust grains are assumed 
to be determined only by a balance of the electron and 
ion currents flowing into the grains. Wang Zheng-Xiong 
et al. studied the effects of ion temperature [14] and the 
negative ions on the zero-point of dust charging (neutral 
point of the dust charge) [15] in the plasma sheath. In 
Ref.[9], Jin. Yuan. Liu et al. studied the distribution and 
suspension of the particles in a plasma sheath by the fluid 
model and self-consistent grain charge variations. In their 
study, only the negative charged particles were taken into 
consideration. 

In this study, we report a numerical method for the so- 
lution of the sheath modiflcation created in the presence 
of cold fluid dust grains with variable charge. In addition, 
we study the dust motion and suspension in the sheath by 
taking into consideration the neutral point of dust charge. 
We shall assume that the charge variation of the dust 
is determined self consistently as in the Ref.[16]. For the 
conditions when the mean free path is long and the dust 
radius is small compared to the Debye length, the orbit 
motion limited theory [17,18] will be used to calculate 
the charging currents. A fluid model and self-consistent 
grain charge variations are chosen to describe the dusty 
plasma sheath, composed of electrons, ions and charged 
dust grains with radius in order of nanometers size. We 
have considered the effect of parameters which modify the 



dusty plasma sheath characteristics, such as the action of 
dust grain radius, the density ratio of ions to electrons and 
the initial velocity of dust. These fundamental parameters 
affect significantly the zero point of the dust charge, this 
neutral point affects effectively the spatial dust density 
distribution. This work is organized in the flowing fashion. 
In sec. 2 The analytical model of the sheath structure and 
the dust charging in the sheath are presented. In sec. 3 
the numerical results are presented and discussed. Finally, 
a conclusion is given in sec 4. 



2. DUSTY PLASMA SHEATH MODELS 

We consider a stationary state [dt = 0), collisionless 
and unmagnetized dusty plasma sheath in contact with 
a planar wall, consisting of electrons, ions and nano 
size charged grains, which has one-dimensional coordinate 
space. Assume the sheath region lies between x = Q (the 
sheath edge), and the wall can be located any where 
in the region a; > 0. At the sheath edge x = 0, the 
potential is assumed to be taken zero, the electron, ion 
and dust densities rieo, n^o, "-do satisfy the quasi-neutrality 
condition e(nio — f^eo) — Qdondo = 0, here e is the electron 
charge and Qdo the dust charge at the sheath edge. The 
electrons are assumed to be in thermal equilibrium state, 
thus the density He is given by the Maxwell-Boltzmann 
distribution. 



Ue = neoexp( — ) 



(1) 



where 4> is the spatial electric potential in the sheath, Tg 
is the electron temperature and n^o is the electron density 
at the sheath edge. 



The ions, arc treated as a cold fluid, governed by the 
continuity and the momentum conservation equations. 



dimensionlcss variables into equations (1) — (9), we obtain 
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where Ui, rn and rm are velocity, density and mass of 
the ions in the sheath, respectively. 

The dust grains are also treated as cold fluid obeys the 
sources free, the continuity and momentum equations are. 



Udnd = UdoUdo 

dud d(j) 

mdUd^T— = -qdi^ 

ox ox 



(4) 
(5) 



where qd = —ezd is the variable-charge of dust grains, 
Zd is the charge number of the dust, md, Ud and Ud are 
respectively the dust particle mass, density and velocity 
of the dust grains. The dust charge arises from plasma 
currents due to the electrons and the ions reaching the 
dust grains surface. In this case, the variable dust charged 
is determined self-consistently by the charge conservation. 



dqd J 
ox 
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(6) 



where le and li are the plasma electron and ion cur- 
rents, respectively, flowing the dust grains surface. The 
charging currents are from the thermal electrons and the 
cold ions [19,20]. 
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3. NUMERICAL RESULTS 

We have numerically solved equations (10) — (16). In 
the flowing calculations, we take some typical plasma 
parameters that are representative of argon plasma, such 
as rieo = 3 * 10^"{cm-^), T^ = 2(eF). The starting point 
is taken at ^ = 0, where we set 77 = 0, -gr = 0.01, 
Vio = 1.1. We consider the dust spherical particles of radius 
is in the range of 50 — 200nTO with uniform mass density 
Pd = 2.2{g/cm^) so that Te/e^ = T^/IA * lO'^ w 1.43 * 
109(m) for Te = 2(ey). 
The results are presented in the following flgures. 
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where r^ is the dust grain radius, me is the electron 
mass, (j)d the dust surface potential. Then, the dust charge 
qd can be determined for a given rd by the dust charge 
surface potential relation (pd 



qd 

rd ' 



The system of equations is completed by the Poisson's 
equation. 



-— = -47r(e(nj - rig) + qdrtd) 



(9) 



In order to solve the system equations (1) — (9), we 
introduce new dimensionless variables 77 = — ^ is the 
normalized sheath potential, the normalized dust charge 

— ^^, the density ratio of ion to electron die = -^^^^^ 
and the density ratio 5de — zr^i as well as normalize the 
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■^ with respect to the electron Debye 
— )i/2_ yJq also normalized ion and dust 



coordinate by ^ = 

length \d^{j^ 

velocities Vd—"^ , Vi — ^, respectively. 

Here Cd = (■^) is the dust-acoustic speec 

Ci — (— ^) is the ion-acoustic speed. Substituting these 
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Fig. 1. The normalized distributions of dust charge (solid 
line) and dust velocity (dashed line) versus ^ (The 
normalized position) in the plasma sheath for param- 
eters Sie = 1.8, Vdo = 18 and ViQ — 1.1. 

Figure 1 shows the normalized dust charge and dust 
velocity distributions versus ^ (normalized position) in the 
sheath for r^ = lOOnm, Sie = 1-8 and Vdo = 18. It is shown 
that the dust charge is negative at the sheath edge, and 
when the dust particles are moving to the wall, the dust 
charge decreases to zero, then close to the wall they are 
charged positively. This is due to the fact that the density 
of electrons falls faster than that of ions in the sheath. 




wall. The curve of the dust charge reveals the existence of 
a neutral point and the dust velocity curve reveals the 
existence of a minimum at this position, this is because 
the dust velocity and dust charge are related in the sheath 
by the momentum equation of dust( Eq.l4). It is clear, 
from this equation, that the dust velocity represents an 
extremum at neutral point of dust charge rjd — 0. 

Figure 2 shows the distribution of dust density in the 
sheath for the conditions same as for fig.l. It takes on the 
oscillatory structures at first, then sharp peak emerges and 
rapidly drops. The result is the same as J.Y.Liu's result[8]. 
The sharp peak emerges at zero point of dust charge where 
the dust velocity is minimum. 



Fig. 2. The normalized dust density distribution versus ^ (normal- 
ized position) in the plasma sheath with the same parameters 
as figure 1. 




Fig. 



3. The normalized dust density distribution versus ^ (normal- 
ized position) in the plasma sheath with the same parameters as 
figure 1 and different values of dust grains radius r^ = 200nm 
(solid line), r^ = lOOnm (dashed line) , r^ = 70nm (dotted 
line) . 



' 


' 


N^, r^=70 nm 






Ni, r^=70rm 




N^, r^=200 nm 


; 
: 




Ni, r^=200 nm 


\ 
t 

i" 




- 


1': ■■. 




- 















Fig. 4. The normalized density distributions of ions and electrons 
versus ^ (normalized position) in the plasma sheath with 
different values of dust grains radius r^ = 200ram r^ = 70nm. 

Then, the grain charge is due only to the ion flux, where 
the electrons vanish and there is a pure ion sheath near the 



Figure 3 shows the dust density distribution under 
various values of the dust grain radius r^, while the 
equilibrium dust velocity is Vdo ~ 18. We notice that, the 
dust density distribution shows a peak, which indicates 
that in this region, more dust particles are gathered. With 
the decrease of dust radius, the peak of the dust density 
moves towards to the wall. This is because the zero point 
of the dust grain charge for large dust grain radius closes 
faster to the sheath edge in opposite to small grain radius. 
Thus, the region of dust disperse becomes wider. 
Figure 4 shows comparisons between density distributions 




Fig. 5. The normalized net space charge density versus ^ (normal- 
ized position) in the plasma sheath with the same parameters as 
figure 1 and different values of dust grains radius r^ = 200nm 
(solid line), r^ = lOOnm (dashed dot line), r^ = 70nm (dashed 
line). 

of electrons and ions for two cases r^ = TOnm and r^ = 
200nm in the sheath. Both the ion density Ni and the 
electron density N^ fall slowly in the case of (r^ = 70nm) 
compared to the case of (r^ = 200nTO). For two cases from 
a certain position the electrons begin to vanish and the 
sheath become an ions. 

Figure 5 shows the net space charge for various values 
of the dust grain radius r^. The net space charge takes the 
oscillatory structures at first, then the positively charged 
gather in certain region, so the net space charge density 
in such region of the sheath emerges at a peak value. The 
case of the dust radius (r^ = 200nm) shows the net space 
charge increases fast at the sheath edge, because the elec- 
tron and negatively charged dust vanish more rapidly than 





Fig. 6. The normalized spatial electric potential versus ^ (normal- 
ized position) in the plasma sheath with the same parameters 
as figure 1 and different values of dust grains radius r^j = 70nm 
(solid line), r^ = lOOnm (dashed lme),r^ = 200nm (dashed dot 

line) . 



Fig. 9. The normalized dust charge distribution versus ^ (normal- 
ized distance) in the plasma sheath with the same parameters 
as figure 1 for various density ratio of ion to electron at the 
sheath edge Si^ = 3.8 (solid line), <5ie = 2.8 (dashed line), 
Sie = 1.8 (dotted line). 





Fig. 7. The normalized dust density distribution versus ^ (normal- 
ized distance) in the plasma sheath with the same parameters 
as figure 1 for various dust velocity at the sheath edge v^q = 20 
(solid line), v^q = 18 (dashed line), v^jq = 22 (dotted line). 




Fig. 8. The normalized dust charge distribution versus ^ (normal- 
ized distance) in the plasma sheath with the same parameters 
as figure 1 for various dust velocity at the sheath edge tJ^o = 20 
(solid line), v^q = 18 (dashed line), v^jq = 22 (dotted line). 



Fig. 10. The normalized dust density distribution versus ^ (normal- 
ized distance) in the plasma sheath with the same parameters 
as figure 1 and for various density ratio of ion to electron at 
the sheath edge Sig = 3.8 (solid line), Sig = 2.8 (dashed line), 
Sif, = 1.8 (dotted line). 

those in the cases of radii (r^ = lOOrirn) and (r^ — 70nm). 
Thus, the net space charge is greater with increasing the 
radius of dust. Accordingly the spatial potential in the 
sheath increases more rapidly for dust of large radii (see 
fig. 6), because of its high net space charge density. 

In figures 7-8, we study the effect of the parameter 
Vdo (the equilibrium dust velocity) on dust charge and 
density distributions. It is shown that the bigger of the 
dust density nearer from the sheath edge. It is also found 
that greater values of Vdo favor the expansion of the region 
of dust disperse, this make distribution of dust approach 
to the wall fig. 7. This is because the neutral point of dust 
charged in the case of Vdo == 18 is closer to the sheath 
edge than that in the cases of Vdo — 20 and Vdo = 22 , 
the reason for it is that the negatively charged particles of 
Vdo = 18 vanish faster than that in the cases of Vdo = 20 
and Vdo = 22 see fig. 8. 



The effect of the density ratio of ion to electron (the 
equihbrium dust density) on the dust density and the dust 
charge distributions in the sheath is studied in the figures 
9-10. By increasing value of (5ie, there are more ions than 
electrons at the sheath edge, so the dust grains are less 
negative and for that the grains charge drops quickly to 
zero. One can see that by increasing the density ratio of 
ion to electron 6ie, i.e much ions at the sheath edge. It 
is shown that for the great value oi Sie = 3.8 compared 
to 6ie = 1.8, the neutral point of dust charge closes faster 
near the sheath edge. Accordingly the position of the peak 
of dust density shifts towards the wall for the small values 
of 6ie fig-9. One can see also that the peak of dust density 
decreases with increasing Sie, this can be explained by the 
fact that at the sheath edge, the dust density n^o is given 
by the quasi neutrality condition (n^o = {Sie — l)/^do), so 
that the normalized density Nd{= nd/rido) is low. 

4. CONCLUSION 

In this paper, we have presented the spatial distribu- 
tion and suspension of the nanodust particles with self- 
consistent charge variation in plasma sheath using fluid 
model. It is found that the presence of dust can lead to a 
very different physical behavior compared to that of the 
electron-ion plasma. The details of the dust particles, such 
as the variations in size, charge and mass can affect the 
spatial distribution of dust. Furthermore, the dust grains 
with different initial velocities, sizes and densities entering 
the sheath, show different motion states: returning at the 
sheath edge and moving at the wall, thus; the region of 
dust disperse is modified. In addition, it is shown that at 
the sheath edge the dust particles are charged negatively 
then, it are become charged positively in the sheath. More- 
over, we also found that the neutral point of dust charge 
shifts by different boundary conditions at the sheath edge. 
The results of the present investigation can be used in 
understanding that the neutral point of the dust charge 
play an important role on the dust particles motion in the 
sheath. 



[15] Wang Zheng-Xiong, Wang Wen -Chun, Liu Yue, Liu Jin- Yuan 
and Wang Xiao-Gang, J. Plasma. Physics 70 (2004) 577-581. 

[16] M. Tribeche, H. Houili, and T. H. Zcrguini, Phys. Plasmas 9 
(2002) 419-429. 

[17] E.G. Whipple, Rep. Prog. Phys AA (1981) 1198-1250. 

[18] J. E. Allen, Physica Scripta 45 (1992) 497-503. 

[19] A. Barkan, N. D'Angelo and R. L. Merlino, Physical Review 
Letters 73 (1994) 3093-3096. 

[20] J. Gorec, Plasma Source Sci. Technol 3 (1994) 400-406. 



REFERENCES 



J. Y. Liu and Y. Liu, J. Apll. Phys 97 



[I] Z. X. Wang, X. Wang 
(2005) 1-5. 

[2] J. X. Ma and M. Y. Yu, Phys. Plasmas 2 (1995) 1343-1345. 
[3] M. Y. Yu, H. Saleema and H. Luo, Phys. Fluids B A (1992) 

3427-3431. 
[4] M. K. Mahanta, K. S. Goswami,Pramana-J. Phys 56 (2001) 

579-584. 
[5] D. Tskhakaya , P. K. Shukla , B. Eliasson, Physics Letters A 

331 (2004) 404-408. 
[6] T. E. Sheridan.J. Appl. Phys 98 (2005) 2-6. 
[7] T. Nitter, Plasma Source Sci. Technol 5 (1996) 93-111. 
[8] G. Arnas, M. Mikikian, G. Bachet, and F. Doveil, J. X. Ma and 

G. X. Yu, Phys. Plasmas 7 (2000) 4418-4422. 
[9] J. Y. Lieu, D. Wang and T. G. Ma, Vacuum 59 (2000) 126-134. 
[10] D. Wang and M. Sheng, J. Appl. Phys 94 (2003) 1368-1373. 

[II] D. Wang and X. Wang.J. Appl. Phys 89 (2001) 3602-3605. 
[12] N.-G. Wang, B.-S. Xie, Phys. Plasmas 9 (2002) 717-720. 

[13] G. Lin, Wci-gui Feng, M. M. Lin, Communications in Nonlin- 
ear Science and Numerical Simulation 13 (2008) 1287-1293. 

[14] Z. X. Wang, J. Liu, Y. Liu and X. Wang, Chin. Phys. Lett, 21 
(2004) 697-699. 



